Chromosome Research

The genus Rhynchospora Vahl (beak-sedges) comprises approximately 381 accepted species with a worldwide distribution, all of which possess holocentric chromosomes, where centromeric activity is distributed almost along the entire chromosome. Despite the recent advances, the mechanisms governing the dynamics of meiotic recombination in holocentric plants remain poorly understood. Here, we developed haplotype-specific oligo-FISH probes for chromosomes 1, 2, and 3 based on a haplotype-phased genome assembly of Rhynchospora breviuscula (n = 5), enabling homolog-specific chromosome painting. Each probe set was labelled with a distinct fluorophore and hybridised in situ to metaphase chromosomes of the reference plant and seven F1 individuals derived from self-crossed reference plants. This approach allowed the unambiguous discrimination of homologous haplotypes and the indirect visualisation of crossover (CO) events in recombined chromosomes. We observed that recombination events were predominantly located in terminal chromosomal regions, consistent across individuals. These results corroborate previous findings from single-cell recombination mapping and provide independent cytological validation of the recombination landscape in this species. Our study establishes haplotype-specific chromosome painting as a robust tool for high-resolution mapping of meiotic recombination in holocentric plants across generations. Furthermore, these probes provided a foundation for future investigations into inverted meiosis, a mechanism characterized by an alternative pattern of chromosome segregation in holocentric species.

Whole-genome sequencing identifies focal DNA amplifications with base-pair resolution but cannot determine whether amplified sequences reside on extrachromosomal DNA (ecDNA, also known as double minutes) or within chromosomally integrated homogeneously staining regions (HSRs). DNA fluorescence in situ hybridization (DNA-FISH) metaphase spreads remain the gold standard for distinguishing these amplification states at single-cell resolution. Here, we present a detailed protocol for DNA-FISH metaphase spreads using human cancer cell lines, encompassing cell culture, metaphase arrest, hypotonic treatment, fixation, chromosome spreading, fluorescent probe hybridization, and fluorescence imaging. The protocol incorporates intermediate quality-control steps to verify successful chromosome dispersion and optimize metaphase spread quality, making the workflow accessible to laboratories without specialized cytogenetics expertise. Results demonstrate clear visualization of ecDNA and HSR amplification states using locus-specific probes and illustrate common technical artifacts that can affect interpretation. This protocol provides a robust and reproducible approach for studying the structural organization of oncogene amplification in cancer cells.

Buzz pollination is an important behavior in which bees use vibrations to extract pollen from poricidal anthers. However, the extent to which vibration frequency influences pollen release remains unclear. Here, we quantified pollen expulsion from Solanum sisymbriifolium anthers subjected to harmonic excitation over a broad frequency range encompassing the anthers first natural frequency. We excited anthers to expel pollen and measured anther kinematics and pollen release using high-speed videography. Particle tracking enabled continuous estimation of pollen release throughout each buzzing event, allowing both initial pollen flux and total pollen released to be quantified. Pollen release depended strongly on excitation frequency. Initial pollen flux, total pollen release, and anther kinematics peaked when excitation frequency approached the anthers natural frequency. Anther tip velocity amplitude exhibited the strongest correlation with total pollen release (r = 0.755) and initial pollen flux (r = 0.898). Experimental observations were compared with nonlinear and linear statistical models of pollen release. While both models captured trends in normalized pollen flux, they overpredicted total pollen release, suggesting that adhesive interactions play important roles during extended buzzing events. These findings demonstrate that anther structural dynamics influence pollen release and suggest that vibration amplification may improve the efficiency of buzz pollination.

Microscopes are essential tools for discoveries on a scale invisible to the unaided human eye. The development of immuno-fluorescence followed by molecular biology techniques and fluorescent fusion proteins have revolutionised the use of optical microscopy in bioscience. The quality of the data produced is dependent upon the sample, its preparation and the instrument used. However, instruments can degrade over time without easily visible changes to the produced images and, in turn, negatively impacts results. By testing instruments and doing comparisons between results over time and between different instruments, problems can be highlighted and corrective action can be taken. Using small fluorescent beads the point spread function (PSF) of the microscope can be recorded and the image resolution measured. Beads were prepared in a concentration matched to the field of view size and dried onto coverslips and mounted on slides. The beads were then imaged as 3D Z-stacks of sufficient size to fully enclose the PSF of the system. This data was uploaded to OMERO and processed using OMERO-metrics, an OMERO plugin developed for this purpose. This paper summarizes the development of workflows and protocols to enable this process, presents the results obtained and demonstrates the detection of significant instrument issues.

Reporters are widely used in Drosophila genetics to visualize gene expression and cell lineages. However, uncharacterized limitations in specific reporter lines can lead to data misinterpretation. Here, we identify a consistent, driver-independent tdTomato signal in the adult proventriculus from the widely used lexAop-tdTomato.nls reporter line. This signal was observed across multiple lexA driver combinations and was directly detectable in lexAop-tdTomato.nls responder-alone adult proventriculi lacking any lexA driver and without antibody staining. In contrast, no comparable native red fluorescence was detected in larval proventriculi under the same no-antibody imaging condition. Mouse and rabbit anti-RFP immunostaining further supported the presence of proventriculus-associated tdTomato/RFP antigen in adult responder-alone animals. In larval responder-alone proventriculi, antibody-amplified staining was antibody-source-dependent: a detectable signal was observed only with rabbit anti-RFP, whereas mouse and rat anti-RFP produced no reliable detectable signal under the same staining condition. A driver-matched comparison using lexAop-RFP.nls did not reproduce the proventricular signal, arguing against detectable ectopic activity of the tested lexA driver in this tissue. However, because lexAop-tdTomato.nls and lexAop-RFP.nls differ in reporter/transgene architecture and possibly genomic insertion context, the underlying cause cannot be assigned specifically to the lexAop sequence. Our findings highlight the necessity of including driver-negative and no-antibody controls when using this reporter line in adult Drosophila proventriculus and gut studies.

In most organisms, all cells inherit the same genome, and many mechanisms exist to preserve its integrity across cell divisions. Programmed DNA elimination (PDE), the targeted removal of specific genomic regions from somatic cell lineages during early embryogenesis, is a striking exception. Since its discovery in parasitic nematodes over a century ago, PDE has been observed in diverse eukaryotes, including ciliates, arthropods, and vertebrates. However, the mechanisms, functions, and evolutionary origins of PDE remain poorly understood. Here, we describe the discovery of PDE in three species of the free-living nematode genus Caenorhabditis. Multiple genomic regions are precisely eliminated from somatic cells during early embryogenesis, resulting in chromosome fragmentation and the loss of key germline genes. The sites of elimination are strongly associated with conserved sequence motifs that likely direct DNA breakage. Comparative analyses indicate that PDE was present in the last common ancestor of Caenorhabditis and subsequently lost early during the evolution of many species, including C. elegans. The presence of PDE in the ancestors of one of biology's most important model organisms, together with recent discoveries in other eukaryotic lineages, reveals PDE to be a far more widespread and significant feature of evolution and development than previously recognised.

High-quality plasmid DNA purification at high throughput remains a significant bottleneck in molecular biology and bioengineering. Current methods frequently fail to deliver sufficient yields of pure, transfection-grade DNA required for genetic engineering applications in mammalian cells. Here, we present a Biofoundry-based automated pipeline using the CyBio FeliX robotic liquid handling platform to rapidly purify plasmid DNA with minimal manual intervention. The protocol leverages Solid Phase Reversible Immobilisation (SPRI)-based magnetic bead technology to ensure consistency, scalability, and DNA purity suitable for downstream viral particle production and mammalian cell transfection. The pipeline supports flexible processing of between 8 and 96 samples per run, making it adaptable across a wide range of experimental scales. The protocol is openly available via Earlham Institute GitHub repository, enabling broad adoption across the bioscientific community and contributing to the growing toolkit of reproducible, scalable engineering biology workflows. In this work, we employed an integrated robotic pipeline to process 528 pooled DNA plasmids and built a Lentiviral DNA plasmid library for lineage tracing, validated the library by sequencing, and demonstrated efficacy in downstream mammalian cell transfection experiments.

Eukaryotic genomes often have fragile sites where chromosomes are particularly prone to break. In Drosophila, when dicentric ring chromosomes try to segregate, they break at nonrandom hotspots. Here, we precisely map breakage hotspots produced by dicentric ring chromosomes in Drosophila. Our study provides three key results about the nature of dicentric chromosome breakage. First, duplications produced by dicentric ring chromosome breakage are surprisingly complex and involve many structural rearrangements, indicating that healing of these breaks is not a simple process. Second, characterization of one particular hotspot showed that new termini all occurred within a single intron of a large testis-expressed gene, suggesting that replication-transcription conflict may be a key determinant of chromosome fragile sites. Third, the new ends are often located near preexisting transposons, suggesting that transposon insertions may contribute to fragility or participate in stabilization of broken ends.

Heterochromatin protein 1 (HP1 is an intrinsic component of heterochromatin domains where it is involved in a diverse set of functions including heterochromatin spreading and organization, chromatin compaction and transcriptional silencing. It has been suggested that HP1 functions through a phase separation mechanism, a process that has been observed in vitro in the presence of N-terminal phosphorylation, nucleic acids and nucleosome arrays. HP1 can also interact with numerous binding partners that contain a specific motif called an HP1 access code (HAC). HACs recognize and bind to an interface formed by the chromoshadow (CSD) domains in the HP1 homodimer, the functional form of the protein. It has been shown that some HP1 binding partners can enhance its phase separation ability while others disrupt the process. Here, we focus on the interactions between HP1 and three binding partners, namely the p150 subunit of the chromatin assembly factor 1 (CAF-1), the N-terminal domain of the lamin B receptor (LBR), and the mitotic protein Shugoshin 1 (Sgo1). Using phase separation assays, we show that CAF-1 prevents HP1 phase separation while LBR and Sgo1 enhance it. Binding assays, mutational studies, NMR spectroscopy and computational analysis allow us to dissect the contributions of the HAC motifs, the charge patterns of the binding partner sequences and the role of N-terminal phosphorylation on HP1 in condensate formation. Our results demonstrate that each binding partner uniquely balances these contributions to modulate the properties of HP1, while electrostatic interactions dominate the regulation of phosphorylated HP1. These results suggest that HP1 binding partners play an important role in the modulation of its properties and the regulation of its functions in distinct biological contexts.

In trypanosmatids genes, transcribed by RNA polymerase II do not have canonical promoters and are organized into directional gene clusters that mature into monocistronic transcripts by a co-transcriptional process known as trans-splicing. Even though gene expression is regulated mainly post-transcriptionally, it is currently understood that chromatin and epigenetics are also involved in this regulation. In eukaryotes, specific signals are normally required for the occurrence of an appropriate transcription initiation. Among them, trimethylation of histone H3 in lysine 4 is the most conserved signal normally detected at transcription start sites of actively transcribed genes. Unlike many model organisms, trypanosomes do not have defined promoters. Instead, transcription initiates in a bidirectional manner from dispersed regions coincident with divergent strand switch regions located between directional gene clusters (DGCs). In T. cruzi, H3K4me3 was observed at the origins of transcription coincident with divergent strand switch regions (dSSRs) in epimastigotes, but it has not been mapped throughout the whole genome at base-pair resolution or in other life stages so far. Here, we set up the CUT&RUN technique for T. cruzi epimastigotes and trypomastigotes. Consistent with a predominant post-transcriptional regulation along the life cycle, we did not find significant differences between life stages. We corroborated that H3K4me3 is enriched at dSSR adjacent to actively expressed DGCs. Moreover, we noticed that this histone mark exhibits different patterns that correlate with the genomic span of the transcription initiation regions and with transcriptional activity. Furthermore, we unveiled that the most actively transcribed DGCs are associated with shorter dSSRs and are located within the core compartment of the genome displaying a more accessible chromatin.

Precise knock-in of fluorescent reporters is a powerful tool for studying the dynamic cellular and molecular processes of embryogenesis. However, conventional CRISPR-Cas9 knock-in of large inserts, such as full-length fluorescent proteins, is inefficient. This has limited its application in many emerging model systems, including sea urchins. Here, we overcome this barrier using a transgenic Lytechinus pictus line that constitutively and ubiquitously expresses a large fragment of mNeonGreen (mNG3K1-10). In this line, fluorescence is only reconstituted when CRISPR-mediated knock-in delivers mNG211, the 11th beta strand of the fluorescent protein, to complement the constitutively expressed fragment. Because this strategy requires integrating only the short 11th-strand, together with short homology arms (~130 nt total), by homology directed repair, it circumvents the size constraints that limit conventional full-length reporter knock-ins using CRISPR. Using this approach, we achieved integration efficiencies of 14-22%, roughly an order of magnitude higher than those obtained with full-length fluorescent protein knock-ins. This provides a streamlined, scalable method for endogenous protein visualization in echinoderm embryos and a valuable resource for studying gene function, morphogenesis, and toxicant response in this classic developmental model.

The haplochromine cichlid fishes of Lake Malawi form one of the most dramatic examples of recent rapid radiation in vertebrates. Here we describe nine new diploid telomere-to-telomere (T2T) genome sequences generated using ultra-long ONT reads, which include 78 ungapped chromosomes. We provide accurate annotations of transposable elements and tandem repeats, identify rDNA cluster regions and putative centromeres, and confirm previously reported large chromosomal inversions. The putative centromeres are primarily composed of satellite tandem arrays of previously reported 237 bp repeats, but notably on most chromosomes these are organised in a novel structure in which four blocks of satellites in alternating orientation are separated by an inverted pair of ~15 kb sequences we term 'centroids', which have similarity to a non-autonomous DNA transposable element and containing potential CENP-B binding boxes. The methylation dip region indicating the likely active centromere always lies between the centroids, whose separation is almost always around 200 kb (interquartile range 151-221kb). A structurally equivalent but non-homologous organisation is seen in the distantly related Etroplus cichlid genera from South Asia. By comparing these structures across chromosomes and species, we suggest how they may have evolved, and potentially how they could contribute to the rampant sympatric speciation seen in these species, based on meiotic drive and chromosome missegregation.

The histone H3K9 methyltransferase Clr4 is essential for heterochromatin formation in Schizosaccharomyces pombe, yet how it searches for and engages chromatin in vivo remains unclear. Using live-cell single-molecule tracking of PAmCherry-Clr4, we quantified how perturbations to Clr4 alter its diffusion, search behavior, and residence times. Chromodomain and SET-domain Clr4 mutants move faster and more isotropically and have reduced residence times at heterochromatin, reflecting impaired substrate recognition. In contrast, deleting Swi6, which has been proposed to sequester Clr4, does not reduce the slow-state fraction or alter Clr4 diffusion, indicating that chromatin engagement is intrinsic to Clr4 rather than HP1-dependent. Anisotropy analysis at short and intermediate displacements indicates that Clr4 does not explore chromatin by simple three-dimensional diffusion but through a guided, distance-dependent search in which it repeatedly samples nearby nucleosomes before disengaging. Across all other perturbations we examined, such as deletion of the CLRC component Rik1, impaired Clr4 ubiquitination, and using cells with an unmethylatable H3K9R substrate, Clr4 dynamics were only modestly affected, and a chromatin-associated population persisted in every background. This robustness indicates that the chromodomain and SET domains are the primary determinants of how Clr4 engages chromatin in vivo, allowing it to continuously sample the genome while maintaining a stable bound population. Our results suggest how the promiscuous sampling of chromatin may also enable Clr4 to establish novel sites of heterochromatin during adaptation.

DNA assembly is a cornerstone of synthetic biology, enabling the construction of bespoke genetic systems for applications ranging from metabolic engineering to DNA nanotechnology. Conventional Gibson Assembly (GA), the most widely used method, relies on 5' exonucleolytic resection and elevated temperatures ([~]50 {degrees}C), which together prevent the retention of 5' modifications and restrict compatibility with temperature-sensitive functionalities. Here, we report a DNA assembly strategy, 3 exonuclease-mediated low-temperature DNA assembly (3LTDA), which generates complementary 5' overhangs while preserving 5' end integrity. This approach enables the efficient assembly of blunt-ended, 5'-functionalised DNA fragments into both linear and circular constructs at ambient temperature (21 {degrees}C), with some assembly observed at temperatures as low as 4{degrees}C. We systematically optimise reaction conditions and demonstrate that this method supports efficient plasmid re-circularisation and multi-fragment assembly, including the construction of a [~]12.5 kbp plasmid from multiple DNA components. Comparative analysis across several DNA substrates shows that, under their respective optimal conditions, this approach matches or exceeds GA performance, improving assembly efficiency by up to 12.8%. Sequence analysis confirms high fidelity with no detectable base-pairing errors across assembled junctions. Crucially, this method preserves chemically functionalised 5' termini, enabling downstream conjugation and biochemical functionality. Retention of azide and biotin modifications was verified through fluorescence imaging, bead-based co-localisation, and enzymatic activity in ELISA-based assays. This is in contrast to GA-assembled controls, which showed complete loss of functionality under comparable conditions. We further assembled 5 kbp dsDNA using 3LTDA from four independent segments, three with different fluorescence reporters, and the fourth containing a biotin group for microparticle conjugation, each on the 5 end. Under fluorescence illumination, bead-bound DNA with all three fluorescence markers were detected. Conventional GA assembled constructs, on the other hand, failed to retain the reporter groups and the fluorescent images did not show the presence of any fluorescent markers. In addition to enhanced performance, the method could also reduce reagent cost and eliminate the need for elevated temperatures, simplifying workflows and expanding the applicability of multi-functionalised DNA constructs. Collectively, this work establishes 3LTDA as a robust, low-temperature alternative to conventional GA, with advantages for applications requiring precise chemical modification, temperature-sensitive components, or deployment outside conventional laboratory environments.

Four-dimensional scanning transmission electron microscopy (4D-STEM) enables the acquisition of diffraction patterns at every probe position in a dense array. For imaging applications this approach offers significant benefits in terms of spatial resolution and contrast enhancement. In this work, we present the development of a synchronous scan generator integrated with SerialEM software to enable automation of complex experimental protocols such as tomography. The proposed hardware functions as an interface between SerialEM, the scan controls of the microscope, a fast annular dark-field detector, and a synchronized trigger for a pixelated detector. Our previous implementation, named SavvyScan, relied on a dedicated computer equipped with a multichannel acquisition and signal-generation cards, as well as a separate microcontroller for synchronization. Here, we report a low-cost implementation based on a Red Pitaya board, utilizing direct programming of its embedded FPGA and Linux server components. We provide detailed instructions for system installation and operation, along with practical guidance for modifying the source code. System performance is validated through oscilloscope measurements and imaging of a replica grating sample. The utility of the approach is further demonstrated by generating a 3D electron tomogram of a cryogenic sample of mitochondria from a tilt series of shadow montage projections.

Mammalian sperm chromatin carries epigenetic information with the potential to influence offspring phenotype, making its faithful characterization essential. It has been suggested that cauda sperm preparations are contaminated by somatic chromatin, that this contamination dominates genome-wide profiles, and that valid results require pretreatment with somatic cell lysis buffer, DNase I, and dithiothreitol. Here we show that properly purified cauda sperm contain no detectable somatic cells or cell-free DNA and that this pretreatment disrupts sperm chromatin organization. SCLB permeabilizes the sperm nucleus, allowing DNase I to fragment the sperm genome in situ, while DTT treatment causes chromatin to leak out of the nucleus. Using ATAC-see, we further demonstrate that Tn5 transposase can access intact protamine-condensed sperm chromatin without DTT, refuting the premise that profiles from untreated sperm reflect contamination. Pretreatment therefore damages the chromatin it claims to purify, and published profiles of untreated cauda sperm are valid and require no systematic re-examination.

Nucleosomes are fundamental units of DNA packaging and gene regulation in eukaryotes. In mammalian sperm, most nucleosomes are replaced by protamines causing extreme chromatin compaction. Various epigenomic studies reported conflicting results on the distribution of residual nucleosomes in mammalian sperm, questioning their potential role in mediating intergenerational inheritance of paternal epigenetic information. Here we performed single-molecule footprinting through Nucleosome Occupancy and Methylome (NOMe) sequencing and applied the Bayesian statistical model nomeR to determine frequencies of nucleosome removal and retention at 103 specific genomic regions in thousands of developing haploid spermatids and mature spermatozoa of mice. While we readily detected footprints of nucleosomes and the transcription factor CTCF in round spermatids, chromatin became transiently highly accessible in elongating spermatids with loss of such footprints, indicating extensive chromatin reprogramming during spermiogenesis. In mature sperm, following nuclear decondensation with recombinant nucleoplasmin, we measured nucleosome occupancy frequencies ranging ~1.2 to 1.7% at mouse loci. In human sperm, nucleosome occupancy varied between ~2.3 to 4.5% at 163 genomic loci profiled. Contrasting mice, chromatin in ~25% of human sperm was accessible upon reducing disulfide bonds between protamines arguing for species specific protamine packaging. Our findings support a stochastic rather than programmed potential role of residual nucleosomes in mammalian sperm in regulating paternal gene expression during ensuing embryonic development.

Multicellularity has evolved multiple times across the eukaryotic tree of life, including among protist lineages. Because transposable elements (TEs) strongly influence genome architecture and gene regulation, understanding their potential impact on genome structure and their relationship with gene expression may provide insight into the evolution of multicellularity. Here, we generated a new genome assembly for the facultatively multicellular amoeba Acrasis kona and performed comparative analyses of TE composition, TE diversity, and TE-density organization across diverse protist lineages. Comparative analyses included unicellular and multicellular representatives from across the tree of eukaryotes, (Heterolobosea, Filasterea, Cristidiscoidea, and Chlorophyceae), including Naegleria spp., Tetramitus jugosus, Capsaspora owczarzaki, Pigoraptor spp., Fonticula alba, Parvularia atlantis, Volvox carteri, and Chlamydomonas reinhardtii. To examine relationships between TEs and gene regulation, we integrated transcriptomic datasets from A. kona, Capsaspora owczarzaki, and Volvox carteri with genome-wide TE-density analyses of differentially expressed genes. TE abundance and composition varied substantially among lineages, with species that exhibit more complex developmental or cellular organization generally containing higher TE proportions than closely related unicellular taxa. Patterns of TE-density organization near up-regulated, down-regulated, and non-differentially expressed genes also differed among systems, ranging from strong TE depletion in A. kona to weaker or cell-type-specific patterns in Capsaspora and Volvox. Together, these findings suggest that transposable elements are associated with multicellularity across diverse protist lineages, although the specific roles they play appear to be complex, lineage-specific, and not yet fully understood.

High-resolution 3D imaging is an important strategy for visualizing and analysing complex skeletal tissue architecture and the bone marrow microenvironment. However, multicolor immunolabeling and imaging of intact skeletal tissues are technologically challenging. The current immunolabeling and clearing methods for intact skeletal elements are very limited, time-consuming and generate low-resolution data or depend on the use of reporter mice. Here, we describe a protocol for efficient clearing and immunolabeling of intact calcified tissues that enables superfast, single-cell resolution, and quantitative 3D light-sheet imaging of intact skeletal elements and teeth. A key aspect of our protocol is the addition of a collagenase digestion step after fixation and decalcification. This step enhances antibody penetration, resulting in deep, comprehensive staining throughout immunostained bones and other calcified tissues. The protocol includes soft tissue removal, fixation, decalcification, bone dehydration, and bleaching, followed by antigen retrieval and permeabilization before the collagenase digestion step. This procedure is performed to prepare the samples for the tissue clearing process that improves bone tissue transparency prior to light-sheet imaging. The entire protocol, from bone collection to image analysis and quantification, takes about 4 days to complete, thus offering significant improvements over previous methods. This protocol is broadly applicable to the visualization of bone microstructure, bone marrow analysis, vascular and neural network mapping, and the study of signaling molecules in bone development and growth. The protocol requires experience with standard tissue processing and immunostaining techniques, and prior experience in tissue clearing and light-sheet imaging is beneficial but not essential. Key pointsO_LIA protocol for efficient clearing and immunolabeling of intact calcified tissues that enables superfast, high-resolution, and quantitative 3D imaging of various intact bones and teeth. C_LIO_LIThe entire protocol takes only 4 days to complete the comprehensive staining and perfect transparency throughout the intact bones, offering significant improvements over previous methods. C_LI Key referencesBiswas, L. et al. Cell 186, 382-397.e24 (2023): https://doi.org/10.1016/j.cell.2022.12.031

Meiosis is a fundamental process in producing both sperm and eggs, yet recombination landscapes often exhibit sexual differences, known as heterochiasmy. Since meiotic proteins are generally expressed in both sexes, the molecular mechanism driving heterochiasmy remains elusive. The -kleisin subunit gene of meiotic cohesin, Rec8, is expressed bisexually in mammals, while its putative teleost ortholog, rec8a, is expressed in a female-biased manner, presumably due to the presence of its paralog originating from the teleost-specific whole-genome duplication (TGD). Here, we elucidated the evolutionary history and expression dynamics of -kleisin genes across teleost lineages. Through comprehensive phylogenetic and synteny analyses, we revealed that major teleost lineages retain two copies of rec8 and rad21, with rec8 loci experiencing drastic chromosomal rearrangements immediately after the TGD. Using in situ hybridization and single-cell transcriptome data in medaka and zebrafish, we demonstrated a conserved sexually biased expression pattern: rec8a is predominantly female-biased, whereas rec8b exhibits male-biased expression during gametogenesis. Furthermore, comparative epigenetic analyses revealed that the conserved sexually biased expression is driven by lineage-specific cis-regulatory elements, rather than conserved ones. Motif analyses imply that regulatory rewiring by transcription factors, including foxl2l in particular, might have played a crucial role in the establishment and maintenance of this paralog divergence. Our findings highlight how whole-genome duplication and subsequent genomic and epigenetic rewiring subdivided the bisexual function of rec8, offering insights into sexually distinct meiotic regulation. HighlightsO_LITeleosts possess a unique -kleisin repertoire originating from the TGD. C_LIO_LITeleost rec8 paralogs exhibit conserved sex-biased expression during meiosis. C_LIO_LIDrastic genomic rearrangements after the duplication rewired the teleost rec8 loci. C_LIO_LIThe conserved expression pattern is governed by lineage-specific CREs. C_LIO_LIThose CREs harbor similar types of TFBSs such as Fox-family TFs. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=94 SRC="FIGDIR/small/731870v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@c8c84dorg.highwire.dtl.DTLVardef@1d65668org.highwire.dtl.DTLVardef@c2d732org.highwire.dtl.DTLVardef@1be54a2_HPS_FORMAT_FIGEXP M_FIG C_FIG